Momentum-space instantons and maximally localized flat-band topological 

Hamiltonians 
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Recently, two-dimensional band insulators with a topologically nontrivial (almost) flat band has 
been studied extensively, which can realize integer and fractional quantum Hall effect in a system 
without an orbital magnetic field. Realizing a topological flat band generally requires longer range 
hoppings in a lattice Hamiltonian. It is natural to ask what is the minimal hopping range required. In 
this paper, we prove that the mean hopping range of the flat-band Hamiltonian with Chern number 
Ci and total number of bands N has a universal lower bound of ^ 4\Ci\/irN . Furthermore, for the 
Hamiltonians that reach this lower bound, the Bloch wavefunctions of the topological flat band are 
instanton solutions of a CP(N — 1) non-linear a model on the Brillouin zone torus, which are elliptic 
functions up to a normalization factor. 



Introduction - The integer and fractional quantum Hall 
(IQH and FQH) effects were discovered in 2D electron gas 
in a strong uniform magnetic field. In the seminal work 
of Thouless et a&, the Hall conductivity quantization in 
integer quantum Hall states was shown to be determined 
by a topological index, the first Chern number C\ of the 
momentum space geometrical gauge field. This observa- 
tion suggests that topological nontrivial bands and quan- 
tum Hall effect can be defined not only in Landau level 
(LL) problems but also in generic band insulators. The 
first band insulator model with a nontrivial band was 
proposed by Haldane 2 . Such a band with nonzero Chern 
number is usually called a Chern band. The systems 
with fully filled Chern bands are called Chern Insulators. 
When a Chern band is flat and partially filled, electron 
correlation effect is important, and topological states 
analogous to FQH states have been proposed^— . Such 
systems are called fractional Chern insulators (FCI), and 
the corresponding FQH effect in lattice models is called 
fractional quantum anomalous Hall (FQAH) effect. 

In these previously studied models, nearly flat bands 
are obtained in tight-binding models with hopping terms 
between nearest neighbor and several further neighbor 
sites^r— . The band width can be tuned to narrower by 
allowing longer range hoppings. To understand the hop- 
ping range condition of a topologically nontrivial flat 
band, in this paper we take an alternative approach: 
we consider Chern insulators with exactly flat nontrivial 
bands, and ask how short-ranged it can possibly be. It is 
conceivable that such a Hamiltonian must contain arbi- 
trarily long range hopping terms to guarantee this exact 
flatness. To characterize how much these Hamiltonians 
resemble the ones with only finite range hoppings, the 
mean hopping range (MHR) of a Hamiltonian is intro- 
duced. We find that the MHR of flat band Hamiltonian 
can be interpreted as the action of a non-linear er model 
defined in the Brillouin zone, and minimizing MHR for a 
given Chern number is exactly equivalent to finding the 
classical solutions of the non-linear a model within the 
given topological sector. The solution of the latter prob- 



lem is known to be the instanton solutions, so that we 
obtain the universal lower bound of MHR and the cor- 
responding flat band Hamiltonians that are maximally 
localized. In the following, we shall first discuss the two- 
band models, and then, generalize the discussion to the 
multi-band case. 

Two-band Models - A generic two-band Hamiltonian 
takes the form: H = cj t /i(k)ck, where Ck has two 
components and /i(k) can be expanded in terms of the 
identity operator 1-2x2 and the three Pauli matrices 
cr x ' v ' z . Assuming that all bands are perfectly flat, we can 
consider ft,(k) of only the following form without losing 
generality: 



h(k) = -n(k)-o-, 



(1) 



where n(k) is a three-component real vector of unit 
length. Clearly, the spectrum is gapped everywhere. 
When the lower band is fully occupied, the system is 
an insulator. The Chern number of the lower band is 
given by 



C\ = — J d 2 k n • (d x n x d y n) , 



(2) 



where d XjV is the derivative with respect to k XlV , namely 
the x and y components of the wavevector k. This for- 
mula implies that the Chern number is equal to the wind- 
ing number of the map n(k) : T 2 — > S 2 , where T 2 is the 
Brillouin zone (BZ). Here and below we will focus on the 
bands with Chern number C\ > since those with C\ < 
can always be mapped to bands with Chern number —C\ 
by taking a complex conjugation h(k) — > h*(k). 

The question we want to address is how localized the 
flat-band Hamiltonian h(k) is in the real space. Thus, 
we write h %j = J -d^ /i(k)e lk '( Rl_R ^, where h+j is the 
coefficient of hopping term from site i located at R^ to 
site j at Rj and V is the volume of the system. Unlike 
the usual cases, the flat-band Hamiltonian does not have 
a hard cut-off in the range of non-zero hopping, hij is 



generically nonzero for any |R$ 



R, 



However, since 



2 



h(k) is gapped and smooth as a function of k, we expect 
hij to be quasi-local in the sense that it decays exponen- 
tially with the increase of | Rj — Rj | . To characterize its 
extent of locality, we introduce the MHR of the Hamil- 
tonian: 

= V ; Tr7,; ; /,„iR, R, ;- ^ k Tr[afe(k)) 2 ] 
Ey^M E k Tr[Mk)^] ' 

(3) 

where the second equality is obtained by Fourier trans- 
form from real space to k-space. For the two-band mod- 
els, we plug Eq. [T]in and obtain 

(R2>= / ||^ n - 9 " n - (4) 

The R.H.S. of the equation takes exactly the same form as 
the action of a 2D 0(3) nonlinear sigma model (NLctM) 
defined in k-space. Also, we notice that the expression 
of Chern number C\ can be interpreted as the instanton 
number in 0(3) NLctM, which provides a lower bound of 
the NLtrM action-^ due to the following identity: 

( R 2)_^ = |0(^n + e^nxd,n) 2 >O (5) 

Thus the minimal MHR is y/ (R 2 ) = v / 2C 1 /tt. The 
lower bound can be reached when the instanton equa- 
tion <9 M n — —e^ v n x <9„n is satisfied. We will refer to 
the Hamiltonians that satisfy these lower bound as max- 
imally localized flat-band Hamiltonians (MLFH). 
The instanton equation is solved by the ansatz: 

n * + m y = TTW' nz = TTW' ( ) 

with <j) an arbitrary holomorphic function of the complex 
variable k = k x + ik y . Due to the periodicity of n(k), 
4>(k) should be a doubly periodic function on the complex 
plane, i.e. cf>(k) — cf>(k + 2ir) = <fi(k + 2m). This type of 
function is called elliptic function in complex analysis. 

The Chern number C\ is simply determined by the 
pole structure of (p(k). As is discussed above, the Chern 
number of the two-band model is the winding number 
of the map n(k) : T 2 — > S 2 , which can be obtained by 
counting the number of times the map sweeps across a 
specific point on S 2 , say the south pole, while taking 
care of the signs of the Jacobian at the preimage of it. 
The south pole of S* 2 corresponds to the poles of <fi(k). 
Thus, the number of poles of <fi(k) in the BZ, counting the 
multiplicity (e.g. an second order pole will be counted as 
2), equals the Chern number C\. Elliptic functions must 
have at least two poles in the BZ, so that the lower bound 
2 1 ix cannot be reached for C\ = 1. 

Multi-band Models - The results above can be directly 
generalized to arbitrary number of bands. A system with 
1 flat bands with energy — 1 and TV — 1 other flat bands 



with energy +1 is described by the following Hamilto- 
nian: 

Mk) = l Nx Ar-2|l,k)(l,k|, (7) 

where |l,k) is an TV-component spinor that denotes the 
Bloch state in the lowest band. Denote \a), a — 1,2,.., TV 
as a local basis in the Hilbcrt space, The Bloch wavefunc- 
tion is denoted by z a {k) = (a|l,k). The TV-dimensional 
normalized spinor z a (k) can be viewed as a CP 1 ^^ 1 
field on the Brillouin zone. The Berry connection is 
a, = — iz^diZ, and the Chern number is defined as 
d = i Jd 2 k e ij diQj. In terms of z(k), the MHR in 
Eq. [3] becomes 

(R 2 > = jL^J d 2 k(D iZ ?(D iZ ), (8) 

with Di — di — idi- Similar to the two-band case, the 
MHR takes the same form as the action of a CP 1 ^^ 1 
model, and thus its lower bound is given by the instanton 
solution: 

<R 2 ) -^^j dMDz)Wz > (9) 

with D = D x + iD y . For N = 2, CP(1) = S 2 and we 
restores the result in the two-band models. 

The MLFH's are determined by the instanton equation 
Dz = d^z — (z^dfrz)z = 0, which is solved by the ansatz: 

z a (k) = f{k,k)<i) a {k), a = l,2,...,JV. (10) 

where cf> a (k) are holomorphic functions and f(k,k) — 
V c f'a(k) < i ) a(k) is a normalization factor (The repeated in- 
dex is summed over). Due to the periodicity of h(k) as 
a function of k, we again find that <^> a (fc)'s are all ellip- 
tic functions. We can always take out an overall factor 
and require (f>\ = 1, so that the independent degrees of 
freedom are N — 1 holomorphic functions, in consistency 
with the N = 2 case in Eq. ©. 

Similar to the two-band case, the Chern number can be 
simply determined by the pole structure of <j> a . As shown 
in Fig. Q] (a), the Chern number can be translated into a 

loop integral C\ — ^ (§qbz + Es fc SQ ) wlt ^ 9BZ 
the boundary of the BZ and C' sa the contours around all 
poles k sa of 4>a- (Without losing generality we assume 
that no pole is at the Brillouin zone boundary.) Due to 
the periodicity of z(k) and thus that of a,(k), the contour 
integral around dBZ vanishes. In the vicinity of k sai we 
have z a ~ (k — k a - l )/\k — k a ' l \ and zp ~ for all other 
components. Evaluating the contour integral around C sa 
we obtain C\ = N a with N a the number of poles of 
4> a , counting multiplicity. 

For C\ < 0, exactly parallel analysis leads to the 
lower bound 4|Ci|/iV7r, and the wavefunctions are anti- 
holomorphic functions of k x + ik y up to normalization. 

Conclusion and discussions - In conclusion, we have 
studied the TV-band models with all bands flat and only 
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FIG. 1. Example of a two-band MLFH with Chern number 
2. (a) The color plot of \4>{k)\ in Eq. ©. The red curves 
stand for the boundary of the Brillouin zone and the contours 
around the poles of <t>(k). (b) The gauge curvature f xy in the 
Brillouin zone. 

the lowest band occupied. We find that the mean hop- 
ping range of such Hamiltonians with Chern number C\ 

has a lower bound ^/(R?) = yfjfa\Ci\- For |d| > 1, 
the bound is saturated by Hamiltonians satisfying an in- 
stanton equation, which have all the components of the 
lowest band Bloch wavefunction as elliptic functions up 
to a normalization factor. The elliptic function form of 



the wavefunctions provides a simple of way of calculating 
Chern numbers, and it may lead to interesting properties 
of the Wannier states and real-space coherent states of 
the Chern bands^^. 

The exact lower bound we obtain is consistent with two 
explicit constructions of almost-flat band models with 
Chern number N proposed recently: the one with N 
bands and only nearest neighbor hopping 14 , and the one 
with two bands and maximal hopping range y/2N—. In 
the continuum limit N — > oo, it is possible to have ex- 
actly flat Chern bands with only local terms, such as 
the Landau level Hamiltonian. However, Landau level 
Hamiltonian is not the N —> oo limit of the MLFH's 
found here, since each Landau level has Chern number 

1. For \C\\ — 1 systems, the MHR lower bound y -^L 

cannot be reached, and it would be interesting to find 
the actual lower bound. We have only considered Hamil- 
tonians with all bands flat. An open question is whether 
smaller hopping range is possible by allowing the unoc- 
cupied bands to be non-flat. 
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